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Abstract 

This  paper  reports  a  modeling  study  of  proton  exchange  membrane  fuel  cell  (PEMFC)  performance  by  using  a  support  vector  machine  (SYM). 
A  PEMFC  is  a  nonlinear,  multi-variable  system  that  is  hard  to  model  by  conventional  methods.  As  regards  the  SVM,  it  has  a  superior  capability  for 
generalization,  and  this  capability  is  independent  on  the  dimensionality  of  the  input  data.  These  two  merits  combine  to  make  it  a  powerful  tool  to 
predict  how  a  PEMFC  will  behave  under  different  operating  conditions.  Herein  a  SYM  model  of  a  PEMFC  system  is  built,  optimized  and  tested. 
First,  the  model  is  determined  with  selected  experimental  data,  and  then  it  is  used  to  predict  PEMFC  performance.  It  is  shown  that  the  model  can 
make  the  prediction  in  10  ms  with  the  squared  correlation  coefficient  as  high  as  99.7%.  Therefore,  the  proposed  black-box  SVM  PEMFC  model 
applies  to  the  simulation,  real-time  control  and  monitoring  of  a  fuel  cell’s  performance. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

A  fuel  cell,  which  can  provide  energy  for  electric  vehi¬ 
cles,  power  plants  and  so  on,  is  probably  the  most  promising 
type  of  renewable  electric  power  source.  The  fuel  cell  vehi¬ 
cle,  for  example,  is  acclaimed  as  the  vehicle  of  the  future  that 
almost  all  manufacturers  have  made  a  developmental  prior¬ 
ity.  Among  various  kinds  of  fuel  cells,  the  proton  exchange 
membrane  fuel  cell  (PEMFC),  with  its  low  operating  temper¬ 
ature,  high  power  density,  high  efficiency,  fast  startup,  quick 
response  and  zero  emission,  is  the  most  suitable  for  vehicles. 
A  sound  model  can  help  size,  simulate,  evaluate  and  optimize 
PEMFC  research  programs  [  1  ] ,  and  is  thus  of  great  importance  to 
researchers. 

The  PEMFC  system  is  a  nonlinear,  multi- variable  system  that 
is  hard  to  model.  To  date,  many  models  have  been  developed,  but 
most  focus  on  the  design  of  the  PEMFC  instead  of  its  application. 
What  matters  most  to  PEMFC  users,  however,  is  not  its  relevant 
internal  details  but  its  performance  under  different  operating 
conditions.  What  they  really  need  is  a  behavioral  model,  with 
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which  they  can  predict  PEMFC  behavior  under  various  operating 
conditions. 

Motivated  by  this  need  and  realizing  the  growing  popularity 
of  the  support  vector  machine  (SYM)  regression  application, 
we  decided  to  make  a  modeling  study  of  the  PEMFC  system 
by  using  a  SYM.  The  support  vector  machine  was  developed 
by  Vapnik  [2]  on  solid  VC-theory  foundations.  It  has  been 
used  for  classification  in  various  domains  of  pattern  recogni¬ 
tion  and  lately,  has  handled  regression  problems  successfully 
[3].  Superior  generalization  performance  is  obtained  from  SVM 
regression  and  more  importantly,  the  performance  does  not 
depend  on  the  dimensionality  of  the  input  data.  In  our  study, 
the  proposed  SVM  model  is  built,  optimized  and  tested  with 
data  obtained  from  a  Ballard  MK5-E  PEMFC  [4].  The  experi¬ 
mental  results  show  that  cell  voltage  is  predicted  with  a  mean 
squared  error  of  0.02%  and  a  squared  correlation  coefficient  of 
99.7%. 

This  paper  consists  of  five  sections.  The  first  section  serves  as 
an  introduction  to  our  study.  The  second  one  discusses  various 
existing  PEMFC  models  briefly.  The  third  one  is  concerned  with 
SYM  theories.  In  the  fourth  section,  the  experimental  process 
of  building,  testing  and  optimizing  our  SVM  PEMFC  model  are 
presented  in  detail.  Conclusions  and  suggestions  for  future  work 
are  summarized  in  the  last  section. 
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Fig.  2.  PEMFC  system. 
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Fig.  1.  PEMFC  components. 


2.  Existing  PEMFC  models 

A  brief  review  of  various  existing  PEMFC  models  is  given  in 
this  section.  Both  advantages  and  disadvantages  of  the  following 
analytical  and  experimental  models  are  discussed. 

A  fuel  cell  is  a  chemical  reactor  that,  when  fed  with  fuel 
and  air,  can  produce  electricity  continuously.  The  fuel  used  is 
determined  by  the  type  of  the  fuel  cell.  In  the  case  of  a  PEMFC, 
pure  hydrogen  is  commonly  used.  Pure  hydrogen  is  fed  into  the 
anode  as  fuel  and  air  is  fed  into  the  cathode  as  oxidant.  Then, 
promoted  by  the  electrocatalyst,  electrochemical  reactions  take 
place,  and  heat,  carbon  dioxide  and  water  are  produced  (Fig.  1). 

Some  analytical  PEMFC  models  have  been  put  for¬ 
ward,  including  simplified  one-dimensional  models  [5,6], 
two-dimensional  models  [7-9],  and  a  more  complex  three- 
dimensional  model  [10].  Although  these  models  can  help  ana¬ 
lyze  and  optimize  fuel  cells  in  the  laboratory  context,  they  have 
limits.  The  simplified  one-dimensional  models  are  short  on  accu¬ 
racy  and  the  two-dimensional  and  three-dimensional  models, 
with  a  few  immeasurable  parameters,  are  too  complicated  for 
most  PEMFC  users.  Consequently,  they  are  not  suitable  for  engi¬ 
neering  applications.  An  empirical  modeling  approach  will  be 
more  practical  for  PEMFC  users,  from  which  they  can  deduce 
PEMFC  stack  responses  without  a  knowledge  of  the  internal 
details. 

There  are  several  experimental  PEMFC  models,  among 
which  the  most  influential  one  was  put  forward  by  Kim  et  al. 
[11] 


C^ceii  =  E0  -  b  log  (  f —  )  -  r-~ - m  exp  («———)  (1) 

\  A  cell  /  A  cen  \  ACell  J 

The  equation  expresses  the  relationship  between  cell  voltage 
cell  and  current  I.  In  it,  Aceii  stands  for  cell  area  and  other 
symbols  are  coefficients:  open  cell  voltage  Eq,  growth  rate  of 
byproducts  n ,  activation  coefficient  b ,  ohmic  resistance  coeffi¬ 
cient  r  and  mass  transport  coefficient  m.  Eq.  (1)  provides  a  good 


indication  of  PEMFC  performance,  but  its  coefficients  relate  too 
closely  to  operating  parameters,  such  as  cell  temperature,  gas 
pressure  and  gas  flow  rate  (Fig.  2).  As  a  result,  a  given  set  of  oper¬ 
ating  parameters  requires  a  corresponding  set  of  coefficients. 
Unless  the  coefficients  change  with  the  operating  parameters, 
the  equation  produces  the  wrong  results.  But  it  is  difficult  to 
keep  the  operating  parameters  unchanged  for  the  whole  period  a 
PEMFC  is  working.  This  drawback  greatly  restricts  the  model’s 
application. 

Jemei  et  al.  [12]  utilized  artificial  neural  network  (ANN) 
methodology  and  came  up  with  an  interesting  and  powerful 
solution,  i.e.  an  efficient  static  PEMFC  model.  In  their  study, 
as  experimental  and  simulated  results  are  very  close,  the  model 
can  be  implemented  without  any  difficulty  in  a  complete  vehicle 
powertrain  simulation.  Highly  efficient  as  it  is,  however,  their 
model  has  two  main  weaknesses.  Firstly,  although  only  four 
experimental  learning  patterns  are  needed  for  the  training  stage 
of  the  network,  the  training  strategy  and  the  topology  of  ANN  are 
normally  determined  by  experience,  which  weakens  the  model’s 
objectivity.  Secondly,  ANN  converges  slowly  and  gets  easily 
trapped  in  local  extremum.  Therefore,  it  is  not  appropriate  for 
real-time  applications. 

In  short,  a  new  modeling  approach  is  needed  to  provide  a  bet¬ 
ter  solution.  In  the  following  sections,  the  S  VM-based  modeling 
method  will  be  presented  in  detail. 

3.  SVM  theories 

The  support  vector  machine,  based  on  statistical  learning 
theories  or  VC-theories,  is  a  novel  and  powerful  tool.  It  was 
originally  developed  at  AT&T  Bell  Faboratories  by  Vapnik  for 
classification  in  various  domains  of  pattern  recognition,  but  has 
recently  expanded  successfully  to  deal  with  regression  prob¬ 
lems.  SVM  has  a  superior  ability  of  generalization  to  ANN.  SVM 
regression  employs  the  structural  risk  minimization  (SRM)  prin¬ 
ciple  to  minimize  risks  rather  than  the  empirical  risk  mini¬ 
mization  (ERM)  principle  that  is  used  in  most  traditional  ANN 
models.  SVM  is  also  more  robust  than  ANN.  In  terms  of  fault- 
tolerance,  ANN  behaves  unsatisfactorily;  all  training  data,  even 
small  noises  and  errors,  can  influence  its  performance.  In  the 
case  of  SVM,  however,  with  the  introduction  of  the  e-insensitive 
loss  function,  it  uses  only  part  of  the  training  data,  i.e.  support 
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Fig.  3.  A  linear  SVM  regression.  Points  lying  outside  the  e-tube  are  named 
support  vectors. 


vectors  (SVs),  to  generate  output.  As  a  result,  the  errors  smaller 
than  s  are  neglected  and  the  notable  errors  are  compensated  by 
slack  variables.  Next,  the  theoretical  background  of  SVM  will 
be  explained  briefly. 

For  a  given  data  set  D  =  {(x/,  yt)\i  =  l,...,/i}cldxR, 
X(  e  Md  is  the  input  data  and  yi  e  R  is  the  output  data.  In  a  linear 
case,  the  mission  of  SVM  regression  is  to  approximate  the  data 
set  D  by  a  function 

/(x)  =  (w,  x)  +  b  with  w g Md,  beM.  (2) 

where  (•,  •)  denotes  the  dot  product  between  w  and  x  in  Md. 
Fig.  3  illustrates  a  linear  SVM  regression.  In  a  nonlinear  case 

ns\ 

(w,  x>  =  -  aJ)K(xi,  x)  (3) 

i= 1 


where  s  >  0  is  a  predefined  constant.  The  regression  algorithm 
uses  the  e-insensitive  loss  function  to  describe  how  the  estimated 
function/(x)  deviates  from  the  true  one.  It  accepts  errors  smaller 
than  s,  but  rejects  larger  ones.  The  regression  can  be  converted 
to  a  convex  optimization  problem 


n 

minimum  <£>(w,  §)  =  \  |  |  w|  |2  +  +  §+) 

i—  1 


subject  to 


Le(y,  /( w,  a)) 


where  C  >  0  is  a  pre- specified  constant,  and  and  are  slack 

variables.  |  |w  1 12  is  the  norm  of  w,  and  minimum  of  \  \  |  w  1 1 2  means 

1  0 

that /(x)  should  be  flat.  With  only  the  constraint  ^  1 1  w|  |  ,  the  opti¬ 
mization  problem  is  sometimes  infeasible;  hence  the  variables 
^  and  are  added  to  promise  the  existence  of  a  solution. 
Besides,  the  addition  of  these  two  variables  allows  of  some 
notable  errors  and  makes  the  algorithm  more  robust.  The  con¬ 
stant  C  adjusts  the  balance  between  the  flatness  and  the  permitted 
errors  beyond  s.  The  area  between /(x)  +  e  and  fix)  —  e  is  called 
the  s-tube. 

Implementing  SVM  algorithm  by  mathematical  program¬ 
ming  is  a  highly  skillful  job,  but  fortunately,  some  software 
packages  have  already  been  developed.  This  paper  employs 
LIBSVM  [13],  an  efficient  and  easy-to-use  software  for  SVM 
classification  and  regression. 


4.  Modeling  PEMFC  by  SVM 


For  a  given  PEMFC  system,  the  relation  between  terminal 
voltage  U  and  current  density  I  is  influenced  by  many  operat¬ 
ing  parameters:  cell  temperature  T,  air  flow  rate  q02,  hydrogen 
flow  rate  glU,  air  pressure  Po2 ,  hydrogen  pressure  Pu2,  relative 
humidity  a ,  membrane  humidity  A,  etc.  Accordingly,  terminal 
voltage  U  of  PEMFC  system  is  given  by 


where  K  (xi,  x)  is  a  predefined  kernel  function,  and  nsw  the  num¬ 
ber  of  support  vectors.  The  values  of  weights  af9  af  are  decided 
during  the  training  process.  The  kernel  function  is  employed  to 
map  the  nonlinear  data  into  a  high  dimensional  feature  space 
where  linear  regression  is  to  be  performed.  Eq.  (3)  is  thus  rewrit¬ 
ten  as 

ns\ 

/(x)  =  yy,+  -  )K(Xi,  x)  +  b  (4) 

i= 1 

Support  vector  regression  uses  the  SRM  principle  to  solve  Eq. 
(4).  The  SRM  principle  requires  that  Eq.  (4)  be  flat,  which  is 
subject  to  the  minimum  of  the  functional 

n 

<P(w,£)=  i||w||2  +  C^-+^+)  (5) 

i= 1 

Constrained  by  the  ^-insensitive  loss  function 


U  =  /(/,  T,  qC>2,  gH2,  Po2,  Fh2?  A,  a,  . . .)  (8) 

Up  to  now,  no  model  has  ever  been  able  to  accommodate  all 
these  operating  parameters.  Our  SVM  model  is  no  exception. 
In  our  experiment,  current  density  /,  which  is  decided  by  the 
uncontrollable  load,  and  cell  temperature  T  are  taken  as  vari¬ 
ables,  and  other  operating  parameters  are  held  constant.  Among 
all  the  operating  parameters,  we  single  out  cell  temperature  T  as 
one  of  our  variables.  The  reasons  for  doing  so  are  as  follows. 

Firstly,  it  figures  significantly  in  determining  terminal  voltage 
U.  Secondly,  it  is  easy  to  measure.  Thirdly,  although  the  stack 
can  be  operated  under  a  constant  pressure  and  at  a  constant 
flow  rate,  temperature  is  hard  to  be  held  constant  in  real-world 
applications.  For  example,  since  the  stack  nominal  temperature 
cannot  be  attained  at  startup  or  when  the  stack  runs  at  a  low 
power  level,  cell  temperature  may  vary  within  a  large  range. 
Finally,  this  simplification  does  not  impair  the  validity  of  our 


0,  if  \y  —  /(x,  o') |  <  s 

| y  —  /(x,  a) |  —e,  otherwise 


Le(y ,  /(x,  Of))  =\y-  /(x,  a)\e  = 


(6) 
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x  Vector  to  be  predicted  x  =  (7, 7) 

Fig.  4.  Illustration  of  the  S  VM  PEMFC  model,  which  is  a  feed-forward  network. 
The  input  is  current  density  and  temperature,  and  the  output  is  voltage.  Support 
vectors  and  weights  are  decided  during  training. 


study,  which  is  aimed  at  modeling  PEMFC  by  SVM.  As 
pointed  out  before,  SVM  generalization  performance  is  inde¬ 
pendent  of  the  dimensionality  of  the  input  data.  This  means 
that  if  our  SVM  model,  obtained  with  the  two  variables,  can 
make  the  prediction  quite  precisely,  then  a  model  taking  into 
account  more  variables  can  do  it  with  the  same  degree  of 
precision.  Besides,  the  new  multi-dimensional  model  can  be 
obtained  easily;  all  needed  is  to  add  more  variables  to  our  model 
and  train  it  again.  Hence  in  our  study  Eq.  (8)  is  simplified 
as 


U  =  /(/,  T)  (9) 

It  follows  that  the  simplified  given  experimental  data  set  is 
D  =  {((/,  T)i ,  Ui)\i  =  1 where  (7,  T)t  e  M2  is  the  two- 
dimensional  input  data  and  77/  e  R  is  the  output  data.  The  aim  of 
our  study  is  thus  to  find  an  SVM  model  that  approximates  Eq. 
(9).  The  scheme  of  our  proposed  model  is  illustrated  in  Fig.  4 
[14]. 

Basically,  it  requires  three  steps  to  build  an  efficient  SVM 
model:  preparing  training  data,  training  the  data  to  obtain  an 
SVM  model,  and  predicting  the  new  input  data  with  the  obtained 
model  [15]. 


Table  1 


An  example  of  cell  temperature  scaled  to  [0,  1] 


Unsealed  (°C) 

Scaled 

20 

0 

24 

0.067 

31 

0.183 

39 

0.317 

56 

0.6 

72 

0.867 

80 

1 

4.2.  Selecting  optimal  SVM  parameters 

The  key  to  obtaining  a  highly  accurate  SVM  estimation  is  to 
choose  a  proper  set  of  meta-parameters  C,  e  and  kernel  param¬ 
eters.  Regression  application  is  essentially  more  difficult  than 
classification.  A  lot  of  work  has  bee  done  on  support  vector 
classification  but  little  on  regression.  Fortunately,  however,  a  few 
recommendations  on  how  to  decide  on  appropriate  SVM  param¬ 
eters,  though  quite  contradictory  and  confusing,  have  been  given 
in  the  literature;  and  several  basic  rules  can  be  extracted  from 
them  [16]. 

The  constant  s  is  used  to  find  the  target  function  that  not  only 
lies  as  close  as  possible  to  the  border  of  the  e-tube  but  also  is 
as  flat  as  possible.  The  larger  e  is,  the  flatter  the  function  will 
be,  and  the  fewer  SVs  will  be.  On  the  other  hand,  however,  a 
larger  s  leads  to  larger  estimation  errors.  Therefore,  the  value  of 
s  ought  to  be  determined  in  a  way  that  it  is  proportional  to  the 
input  noise  level  a. 

The  constant  C  determines  the  balance  between  the  com¬ 
plexity  or  flatness  of  the  function  and  the  amount  up  to  which 
deviations  larger  than  e  are  tolerated.  Because  both  C  and  e  can 
affect  the  complexity  of  the  function,  they  should  be  adjusted  at 
the  same  time.  This  increases  the  difficulty  of  selecting  the  most 
appropriate  parameters.  The  principle  for  deciding  the  value  of 
C  is  that  it  is  equal  to  the  output  range. 

It  is  also  important  to  choose  an  appropriate  kernel  function 
and  then  assign  proper  values  to  its  parameters.  There  are  many 
types  of  kernels:  polynomial  kernel,  hyperbolic  tangent  kernel, 


4.1.  Preparing  training  data 

Training  data  ought  to  cover  the  entire  expected  range 
of  operations  of  the  final  SVM  model.  For  example,  tem¬ 
peratures  from  20  to  80  °C  and  current  densities  from  0  to 
800mA/cm2  all  should  be  covered.  In  most  cases,  training 
data  should  be  scaled,  normally  linearly,  to  [0,  1]  or  [— 1, 
+1].  An  example  of  scaled  temperature  is  shown  in  Table  1. 
Scaling  can  increase  the  training  speed  and  assist  in  select¬ 
ing  optimal  SVM  parameters.  Another  matter  worthy  of  note 
is  that  training  data  and  testing  data  must  be  scaled  to  the 
same  interval.  In  this  paper,  U  —  I  data  obtained  at  24,  31,  39 
and  72  °C  are  used  as  training  data,  and  all  those  obtained 
at  56  °C  as  testing  data  (Fig.  5).  All  the  data,  including  aver¬ 
aged  cell  voltage,  current  density  and  temperature,  are  scaled  to 
[0,  1]. 


Fig.  5.  Characteristics  of  a  Ballard  5kW  MK5-E  PEMFC.  It  is  composed  of 
36  cells;  each  cell  has  a  232  cm2  active  area,  graphite  electrodes,  and  a  Dow 
membrane.  Air  pressure  Pq2  and  hydrogen  pressure  Pho  are  both  regulated  to 
3  atm  [4] . 
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Radial  Base  Gaussian  Function  (RBGF)  kernel,  B-splines  ker¬ 
nel,  two-layer  neural  networks  kernel,  etc.  Polynomial  kernels 
are  given  by 

K(x,  Xi)  =  ((x,  xf)  +  c)p,  pe  N,  c>0  (10) 


Hyperbolic  tangent  kernels  are  given  by 


K(x ,  Xi)  =  tanh(0  +  0(x,  Xi)) 
RBGF  kernels  are  given  by 

K(x,  X[)  =  exp 


(11) 


(12) 


Among  all  the  kernels,  RBGF  kernel  function  is  the  most  popular 
one  that  ought  to  be  tried  first  [17],  and  is  thus  used  in  our  study. 

To  determine  the  optimal  parameters,  we  employ  the  cross 
validation  method,  which  is  often  used  in  practical  applications. 
The  method  consists  of  three  steps. 


(1)  Using  four  groups  of  data  obtained  respectively  at  24,  31, 
39  and  72  °C  as  the  training  data  set. 

(2)  Training  SVM  with  some  parameters  and  predicting  other 
data,  and  then  obtaining  the  degree  of  prediction  accuracy. 

(3)  Changing  the  parameters  and  repeating  the  above  two  steps 
until  a  high  degree  of  accuracy  is  obtained. 


With  RBGF  kernel,  s  is  set  to  0.005.  After  the  above  three 
steps  are  carried  out,  the  final  optimal  C  is  found  to  be  15  and  a 
in  the  kernel  is  0.5. 


Fig.  7.  Averaged  cell  voltage  predicted  by  the  SVM  model  from  20  to  80  °C. 
Current  density  varies  from  0  to  700mA/cm2.  All  the  data  have  been  scaled  to 
[0,  1]. 


4.3.  Predicting  with  the  SVM  model 

After  training,  an  SVM  model  is  obtained,  which  can  be  used 
to  predict  new  input  data.  First,  the  cell  voltage  at  56  °C  with  the 
current  density  in  the  range  from  0  to  700  mA/cm 2  is  predicted. 
A  comparison  between  the  predicted  data  and  the  experimental 
data  is  then  made  to  evaluate  the  model’s  prediction  precision 
(Fig.  6).  Next,  all  the  data  in  the  temperature  range  from  20  to 
80  °C  and  the  current  density  range  from  0  to  700  mA/cm2  are 
estimated  (Fig.  7). 


Fig.  6.  Averaged  cell  voltage  predicted  by  the  SVM  model  at  56  °C.  Current 
density  varies  from  0  to  700 mA/cm2.  All  the  data  have  been  scaled  to  [0,  1]. 
The  mean  squared  error  is  0.097%  and  the  squared  correlation  coefficient  is 
98.6%. 


Fig.  8.  Averaged  cell  voltage  predicted  by  the  new  SVM  model  at  56  °C.  Current 
density  varies  form  0  to  700  mA/cm2 .  All  the  data  have  been  scaled  to  [0, 1] .  The 
mean  squared  error  is  0.02%  and  the  squared  correlation  coefficient  is  99.7%. 

Although  the  results  are  quite  satisfactory,  the  SVM  model 
can  still  be  improved.  As  shown  in  Figs.  6  and  7,  there  are  notable 
errors  between  the  actual  voltage  and  the  predicted  voltage  at 
the  beginning  of  the  curve.  The  reason  is  that  while  the  curve 
changes  greatly  at  the  beginning,  the  training  data  are  too  few  to 
represent  the  trend  of  this  sharp  turn.  These  few  data,  however, 
are  regarded  as  exceptional  data  points  by  SVM  algorithm.  For 
example,  the  first  point  in  Fig.  6  is  treated  as  an  error  and  is  com¬ 
pensated  by  slack  variables.  To  describe  the  curve  in  accurate 
detail,  we  need  more  data  points  for  the  beginning  of  the  curve. 
Therefore,  several  data  points  are  supplemented  to  the  beginning 
of  the  curve  and  the  new  SVM  model  is  trained  again.  Then  the 
cell  voltage  at  56  °C  is  predicted  with  the  new  model  (Fig.  8). 
Comparing  Fig.  6  with  Fig.  8,  we  can  see  clearly  that  the  preci¬ 
sion  at  the  beginning  of  the  curve  is  greatly  improved. 

5.  Conclusions  and  suggestions  for  future  work 

An  off-line  modeling  study  of  a  fuel  cell  using  a  SVM,  which 
takes  PEMFC  as  the  example,  is  reported  in  this  paper.  It  is  shown 
that  the  optimized  SVM  model  is  more  attractive  and  more  com¬ 
petitive  than  other  modeling  solutions  in  that  it  possesses  a  high 
degree  of  precision  and  does  not  require  a  pre-knowledge  of  the 
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fuel  cell.  In  our  study,  training  a  S  VM  with  optimized  parameters 
needs  1583  iterations  and  only  30  ms  on  a  PHI  800  MHz  com¬ 
puter.  The  time  for  prediction  is  no  more  than  10  ms.  Hence, 
this  black-box  SVM  model  is  most  applicable  for  simulation, 
real-time  control  and  monitoring  of  fuel  cell  performance. 

Among  all  the  operating  parameters  that  can  affect  the 
PEMFC  performance,  only  current  density  and  temperature  are 
included  in  our  model.  This  is  because,  in  accordance  with  SVM 
theories,  we  can  incorporate  any  other  operating  parameter  into 
a  more  complete  model  without  degrading  the  model’s  gener¬ 
alization  performance  or  complicating  the  modeling  process. 
Therefore,  an  online  training  and  prediction  algorithm  with  more 
variables  can  be  considered  in  future  work. 
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